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Abstract 

The evolution of two-mode Gaussian state under symmetric amplification, non-symmetric damping and 
thermal noise is studied. The time dependent solution of the state characteristic function is obtained. The 
^s^j ■ separability criterions are given for the final state of weak amplification as well as strong amplification. 
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1 Introduction 

In all practical instances the information and entanglement contained in a given quantum state of the system, so 
precious for the realization of any specific task, is constantly threatened by the unavoidable interaction with the 
environment. Such an interaction entangles the system of interest with the environment, causing any amount 
of information to be scattered and lost in the environment. The overall process, corresponding to a non unitary 
evolution of the system, is commonly referred to as decoherence. To overcome the loss, parameter amplifier is 
added to the system. We in this paper will treat the simultaneous actions of amplitude damping and parameter 
, amplification to two-mode Gaussian state. 

The density matrix obeys the following master equation pQ -£ = — i[H, p] + £p,with the quadratic Hamil- 
tonian H — fo^2j k | (rjjkOjOf. — r l*jk a j a k)i where rj is a complex symmetric matrix (parameter amplifier ma- 
Oh. trix). L\p = ~2~{(nj + l)L[aj]p + njL[a,j]p}, where the Lindblad super-operators are defined as L[o]p = 
— cftop — po^o , Tj is the amplitude damping coefficient of jth mode, rij is the average thermal photon 
number of the environment. Any quantum state can be equivalently specified by its characteristic function. 
Every operator A G B(H) is completely determined by its characteristic function xa '■— tr[AT>(p)\ |2j, where 
T>{p) — exp(/ia^ — p*a) is the displacement operator, with p — [pi, p2, • • • , Ps] , a — [ai, Oa, • ■ ■ , a s ] T and the 
total number of modes is s. It follows that A may be written in terms of \A as j^J: A — f[Yli ^ ± }xa(p)'^ > {—p)- 
The density matrix p can be expressed with its characteristic function X- X = tr[pT)(py\ . The master equation 
can be transformed to the diffusion equation of the characteristic function, it is 0] [5] 
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dt Y d"* d rt 2 < 9 M 



2 The parametric amplifier and the amplitude damping 

The solution of the diffusion equation of the characteristic function can be completely worked out for Gaus- 
sian state in the case of real parameter amplifier matrix rj. We will consider real 77 in the following. If the 
initial state is Gaussian, its characteristic function has the form of x(m,A < *,0) = exp[^m^(0) — p*m T (0) — 
\{p, — p*)"f(0)(p*, — p) T ], the state will keep to be a Gaussian state in later evolution, where m is the first 
moment and is irrelevant to entanglement, 7 is the complex correlation matrix (CM). The time evolution of the 
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complex CM for real amplifier matrix r\ is jSJ 



7(*) = 



M -N 
-N M 



(7(0) 



a (3* 
(3 a* 



M -N 
-N M 



a (3* 
(3 a* 



(2) 



rjN-^M,^ = - V M-^N,withT 



where M and N are the solutions of the following matrix equations 

diag{Ti,T 2 , ■ ■ ■ ,T S }. The solution is M = ±[exp(-r)t-^) + exp(r)t N = |[exp(-7}t- ^) -exp(r)t- S)]. 
The constant matrices a and (3 in have the behaviors — a, (3 = /3 T , they are the solutions of the following 
matrix equations 



2(rja + a*r)) -T/3 - /3T = 0, 
Ta + «r -2^/3- 2/3*77 -r(n+i)-(n + i)r = 0. 



(3) 
(4) 



where n = diag{rii,n 2 , ■ ■ ■ ,n s }. The one mode solution has been known for a long time (see ..6^ and references 
therein) . 

For the two-mode situation, the real amplifier matrix 77 = 77000 + r\\<7\ + 77303, where <tq = I 2 , cti, 03 are Pauli 
matrices. M and N can be simplified to 



M = 



-e~ Clt [cosh(5i/j)cr - sinh( J Bi7j)"CT • fa] + -e^ C2t [cos\\{B 2 t)<j - smh{B 2 t)~a ■ fa], 



(5) 



N = ie~ Cl *[cosh(Si/j)CT - sinh(B 1 t)~a ■ fa] - ^ C2t [cosh{B 2 t)<j a - svah{B 2 t)~a ■ fa]. (6) 



where C h2 = ±m + |(Ti +r 2 ); B 1<2 = yjrjf + (K r i - r 2) ± %) 2 ; b 1,2 = {±r ]1 ,0,±i l3 + \{T 1 -Y 2 ))/B lt2 . We 
will consider the case of symmetric noise, that is, n = nol2.The solutions of equations © and l@J are given in 
the appendix. 



3 The inter-mode amplifier 

The algebra equation of a and [3 in two mode system is complicated in general situation. To investigate the 
entanglement property of the amplifier, we will first consider the case of 770 = 773 = which corresponds to 
inter- mode amplification alone. Thus 77 = 771 a. The solution is (see Appendix) 

01 = i j-77f [a-rfVo + r^Va], (7) 

P - 1 r /2 _/2 ( 8 ) 
1 1 3 71 

where ?7 = no + |, T 3 = rs/To,^ = 2?7i/ro. In the assumption of 7/0 ="3 = 0, we have C\ = C 2 = To/2; 
B X =B 2 = J (r 3 /2) 2 + 7/2 = fcr /2, with fc = ^/V' 2 +r)[ 2 . Denote i' = r t/2,hence 

r' 

M = e _t [cosh(fei')°"o - sinh(fct')-^cr 3 ] (9) 



AT = -e- t 'sinh(A:t / )v-cri (10) 

For the case of weak amplifier, k < 1, that is , rj\ < j(Tq — r 2 ), when i — > 00, we have M,N — > 0. The 

" a /?* " 



state will tend to a Gaussian state which is characterized by the residue complex CM 7(00) 
Peres-Horodecki criterion for separability [E] will be jSJ 



/3 a* 



The 



det7 a det7b + (i - |det7 c |) 2 - tr(7 a o-37 c <T37 b cr37 ( tcr3) > i(det7 a + det7 & ), (11) 
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Figure 1: The separability border of the amplification and non-symmetric damping Gaussian system at t — > oo. 
The weaker amplification range is defined by r§ + A-q 2 < T§. 



where 



a a a c 
a* a b 



, P = 



Pa Pc 
Pc Pb 



Oil P* 

Pi a* 



i = a.b, c. 



(12) 



Then a a = ^(1 - If + I>f ), a b = ^(1 - Tf - T' aV [ 2 ), ^ = 0; p c = j^ V [ (l - if) , p a = fa = 0. The 
state is a x-p symmetric Gaussian state whose Gaussian relative entropy of entanglement can be obtained 
[TU| . The separability criterion now is a 2 a a b + (| — /3 2 ) 2 — 2a a a6/3 2 > |(a 2 + a 2 ), which can be reduced to 
(«a - §)(«& - |) - Pc > 0: tnat is 



[i-r: 3 2 (2no + i) 2 K 2 <4^(i-r^). 



(13) 



For the case of strong amplifier, k > 1, that is , 77^ > jj(r 2 — r|), suppose the complex CM is j(t) 
a' &* 1 

,„ at time i, a direct calculation shows that a = diag{a' a , a b }, p = P c <?\ with 



a 'b 

P'c 



a« + M?(~ - "a) + iV 2 (i - a&) + 2M a N c p c , 
a b + M b 2 (i - otb) + iV 2 (i - a.) + 2M b N c p c , 
p c - M a N c (\ - a a ) - M b N c {\ - a b ) - (M a M b + N 2 c )p c , 



(14) 
(15) 
(16) 



where we have denoted M — diag{M a , M b }, N = N c a\ \ the vacuum initial state is assumed. The state is still a 
x-p symmetric Gaussian state. The separability criterion is (a' a — h)(&'b ~ 5) — P'c ^ 0: which can be written as 

tf{{K 2 1){K 2 - 1) - {K,K 2 - l) 2 (2n + l) 2 r 3 2 } 

-r?f {{k x k 2 i) 2 (i - r 3 2 )r 3 2 (4n 2 - 1) + An 2 (K 2 - l)(Kl 1) 
-4n r?[x 2 - r 3 2 - 2^X2(1 - r 3 2 ) + k 2 (i - k 2 t' 2 )}} - m 2 (K 2 - \){k 2 2 - i)(i - r 3 2 )r 3 2 



> 



(17) 



where K\ — e^ k 1 - )t , K% — e ( fe+1 ) t . When t' — > 00, the separability criterion will be 

V? < 2«o (no + T 3 2 + ^n 2 + (2n + l)r 3 2 ) . (18) 

Inequalities (|13(l 1)18(1 are displayed in Fig. 1 in a combined form. The critical noise no is shown as a function of 
t( x and r 3 . 
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Figure 2: The separability border of the symmetric single-mode amplification, inter-mode amplification and 
non-symmetric damping Gaussian system at t — > oo. The weaker amplification case. 

4 The symmetric amplifier 

The system may undergo symmetric single mode amplification as well as the inter-mode amplification, that 
is 773 = 0, 170 7^ 0. We consider the situation of weak amplification, that is, Ci > B\ (rj > 0, T 3 > is 
assumed). When t — > 00, we have M — > 0, N — > 0, the final state is specified by the residue matrices a and (3 
(see Appendix), ft seems that the separability criterion might be very complicate, however, a direct calculation 
shows that the condition can be written as a quadrature form of the square of the inter-mode amplification 
parameter 7/1, 

S2V1 + sivl + s o > 0, (19) 

with s = (1 - r? 2 ) 2 [r? 4 + 8(1 + If K 2 n (n + 1) + 16(1 - r 3 2 ) 2 n 2 (n + l) 2 ], s 2 = [1 - r? 2 - r 3 2 (2n + l) 2 ] 2 , 
si = -2^ 6 - 8r? 4 n (n + 1) - 2r 3 2 77 4 (8n 2 + 8n + 1) -4(1 - lf)(l - r 3 2 (2n + l) 2 )(2n 2 + 2n + 1) +2^ 2 [8n 2 + 
8n + 3 - 4r 3 4 n (no + l)(2n + l) 2 +r 3 2 (16n* + 32n^ + 24n 2 + 8n - 1)]. 

The border of the separable state set and entangled state set is shown in Fig. 2 with 770 = 0.5 , where only the 
case of weak amplification is shown. Our numerical result shows that the range (in terms of relative asymmetric 
damping quantity T 3 = T 3 /ro = (Ti — r2)/(ri + T2) and the noise no) of weak amplification shrinks as 770 
increasing, the weak amplification entanglement can only be possible when the noise is less than 1/2 photon 
number. 

5 Conclusion 

We have studied the evolution of two-mode Gaussian state under non-symmetric damping, symmetric amplifi- 
cation and thermal noise. The non-symmetric damping is the most general damping of two-mode system. The 
amplification is limited the symmetric case for simplicity, although the most general case of 770 7^ 0, 771 ^ 0, 
773 ^ is also solvable. The case of inter-mode amplification alone is especially simple, its separability criterions 
of final states in both weak and strong amplifications were obtained. The separability criterion of the final state 
of symmetric amplification (770 ^ 0, 771 ^ 0, 77 3 = 0) is given for weak amplification. Inter-mode amplification 
parameter 771 is crucial for entanglement. 

In the weak amplification case, final state entanglement is only possible when the thermal noise uq is less 
than 1/2 photon number. When the single mode amplification parameter 770 increases, the entanglement range 
in terms of relative asymmetric damping quantity T' 3 — r 3 /ro = (Ti — iy/^i + T2), the noise no and inter- 
mode amplification normalized parameter r][ shrinks. Lower T' 3 , no and higher 77^ are required for the state to 
be entangled when 770 increases. 
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Appendix: The residue matrices a and (3 

In the two mode situation, denote a = Yli=o a i a ii au a i are rea l due to = a; denote = J2i=o 1 3 A 17 *' Pi = 
0iR + i0u, the C2 item is nullified due to T = 0. Let r ,3 = ^ (Fi ± T 2 ), then T — T ao + r 3 (T3. Together with 
j] = Y]j—q i 3 Vi^i an d n = ^oco, a ll the matrices in Eqs.JSJ) (@J) are expressed in the basis of Pauli matrices. By 
comparing the coefficient of the Pauli matrices, from Eqs.© I@J, we obtain two groups of equations, the first 
group equations containing a =(a , ai, a 3 ) T , R = (0or, 0ir, 03r) T are 



with 



Ga- 


■E0r 






-i)(r ,o,r 3 ) T 


Ea- 


-G0 R 













' To 





r 3 " 






Vo Vi V3 
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r 





, E 
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Vi Vo 




r 3 





To 






V3 f]Q 



The second group equations containing (012, Poi, 0n, 03i) have a solution (012, 0oi, 0u, 03i) = 0. The solution 
to Eqs.Q I© is 

a = (fi + ^)(G-EG- 1 E)- 1 (T ,0,T 3 ) T 
= G~ l Ea 

When 773 = 0, the solution is 

+4 Vom [(2T 2 - r 2 )(r 2 - r 2 ) + 4r|(r 7 2 - v 2 ) - srfa 2 }^ 
+r r 3 [i6(2 % 2 - rii){ n l - v 2 ) + a v 2 (y 2 - rg) - sr 2 ^}, 



= (no + i)A- 1 {2r ry (r^-4r /o 2 )[r2-r| + 4(^- % 2 )]a 

+2T om [(T 2 - T 2 ) 2 + 8 V 2 (2 V 2 - 2 V 2 - T 2 ) + 4 V 2 (r 2 - T 2 )}^ 
+27 ?0 r 3 [16(r 7 2 - rjl) 2 + T 2 (r 2 - r 2 - 8?y 2 ) + 4r 2 (?7 2 - rft)]a 3 }, 

where A = (r 2 - 4 % 2 )[(r 2 — T| — 4 % 2 - A V 2 ) 2 - 4ry 2 (r 2 + 4ry 2 )] = (r 2 - 4r ?0 2 )[(r + 2,7 ) 2 - (r§ + 4ry 2 )][(r„ 
2^o) 2 — (r 2 + 4?7 2 )]. When 770 = ^3 = 0, the solution is 

• - ( " + ^ [r (r 2 -r 2 )a + 4r 3?7 2 a3], 



r Q (r 2 -r 2 -4 V 2 ) 
2(n + i)r ?1 (r 2 ~r 2 ) 
p r a (T 2 -T 2 -A V 2 ) ai - 
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